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Abstract 

' We reconsider the problem of modelling static spherically symmet- 

ric perfect fluid configurations with an equation of state from a point 
f"^ , of view of that requires the use of the concept of principal transform of 

I a 3-dimensional Riemannian metric. We discuss from this new point 

■ of view the meaning of those familiar quantities that we call density, 

. pressure and geometry in a relativistic context. This is not simple 

\ semantics. To prove it we apply the new ideas to recalculate the max- 

O^' imum mass that a massive neutron core can have. This limit is found 

^H^. to be of the order of 3.8 Mq substantially larger than the Oppenheimer 

and Volkoff limit. 
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Introduction 

We review in Section 1 the basic equations of the models being considered 
as well as the concept of principal transform of a 3-dimensional Riemannian 
metric which is at the core of our new point of view to understand these 
models. 

Section 2 is devoted to lay down the fundamental system of equations to 
be integrated. We use spherical space coordinates of the quo-harmonic class 
which allow to implement class smoothness across a sharp border when 
there is one. 

In Section 3 we define the concept of proper mass Mp to be used to define 
the binding energy of the models as the difference between the active 
gravitational mass Ma and Mp. We define also the concept of proper mass 
density which is a fundamental hybrid concept related to Mp and the principal 
transform of the quotient space metric. Its relation with the pressure p 
characterizes the fluid source independently of the solution being considered. 
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Section 4 is devoted to establish that the binding energy Et, can be ob- 
tained as the integral over all space of a localized energy density a which 
depends only on some of the gravitational potentials and its first radial 
derivatives. 

In Section 5 we linearize the fundamental system of equations to obtain 
the linearized expression of the binding energy Eb, which coincides with the 
familiar Newtonian one, thus providing a partial justification to the definition 
of Eh in the non linear regime. We obtain also the linearized expression of 
the energy density cr. 

The last section contains our proposed application of the new point of 
view to the study of massive neutron cores. The equation of state is the usual 
one for a degenerate neutron gas, but both the density and the pressure are 
variants of those used by Oppenheimer and Volkoff [|l[] . Our main result is 
that the maximum limit mass that a neutron core can have is approximately 
3.8 Mq instead of 0.7 Mq in 0. 

1 Static Spherically Symmetric Models 

We shall be interested in this paper on global spherically symmetric models, 
which we shall write using a time adapted coordinate and polar-like space 
coordinates: 

ds^ = -A^r^dt^ + df (1) 

where: 

ds^ = B\r)dr^ + B{r)C{rydQ'^, d^^ = dO^ + sin^ edip\ (2) 
solution of Einstein's field equations: 

Sal3 = -XTaf3, X = 8ttG / (3) 

where the r-h-s describes a compact fluid source, or with fast decreasing den- 
sity, with two "flavours": isotropic or a special kind of anisotropic pressure 
to be presented in a moment. 

The quotient 3-dimensional metric (Q) can be written using a variety of 
supplementary coordinate conditions belonging to two different types: alge- 
braic or differential. The most often used is the curvature condition which 
uses a radial f coordinate with the algebraic condition: 

C = B-^ (4) 
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We shall refer also to any other quantity which assumes the use of this radial 
coordinate with a tilde overhead. 

We shall use here almost exclusively the quo-harmonic condition which 
restricts the r coordinate with the differential condition: 

C' = -{B-C) (5) 
r 

where the prime means derivative with respect to r. The use of a differential 
condition makes possible the construction of global models with a sharp 
boundary with vacuum, something which is not possible when using the 
coordinate f. 

Notice that if B{r) and C{r) are known the curvature coordinate r is 
simply the following function of r: 

f = ^B{r)C{r)r (6) 

while, on the contrary, to obtain explicitly the inverse function is in most 
cases of interest impossible or very cumbersome. 

Despite the emphasis we put on the use of the quo-harmonic coordinate r 
let us be clear from the beginning that the main conclusion of this paper will 
not owe anything to a particular choice of coordinates. It will owe instead 
all to a new concept: that of a principal transformation of a 3-dimensional 
Riemannian metric such as (^. 

By definition, in the particular case we are considering, the principal 
transform of (0) is a new 3-dimensional Riemannian metric Q: 

ds^ = ^^{r)B^{r)dr^ + '^^{r)B{r)C{r)r^dn^ (7) 

such that: i) 

i?;.,; = 0, ^,j,A;,/ = 1,2,3 (8) 

and: ii) 

(q, - T],)g^^ = 0, (9) 

where, with otherwise obvious notations, the quantities with a hat overhead 
refer to the metric (0) and the quantities with a bar overhead refer to the met- 
ric (^. Notice that both conditions above being tensor conditions under any 
transformation of space coordinates the concept of Principal transformation 
is intrinsic to the Killing time congruence we are considering. 

-^See for instance references Q or Q 
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One of the practical conveniences of using polar quo-harmonic space co- 
ordinates from the outset in the very process of model-building is that, when 
appropriate boundary conditions are taken into account, the principal trans- 
form of (0) is just: 

ds^ = dr^ + r^dn'' (10) 

and: 



$(r) = l/B(r), ^(r) = l/^B{r)C{r) (11) 

We shall consider two types of energy-momentum tensors: i) those de- 
scribing standard perfect fluids with isotropic pressure: 

Too = Tij = pQij, (x° = ct) (12) 

and a new type of fluid where the isotropy of the pressure is meant in the 
sense of the principal transform (0): 

Too = Tij = pgij (13) 

i.e. in the sense of ( p!0D if polar quo-harmonic coordinates are used. 

Notice that this second case, that we shall present in Section 3 as being 
the truly isotropic fluid, can be considered from the standard point of view 
as an anisotropic fluid of a particular type. Namely one for which the radial 
and tangential pressures are related to a single function p{r) as follows: 

Pr=pB-\ pt=p{BC)-' (14) 

We can deal with both cases at once introducing in the fleld equations a 
two- valued flag / = or / = 1 and the following expressions for pr and pf. 

p, = (fB-' + l-f)p, p, = (f(BC)-' + l-f)p (15) 

thus / = corresponds to the standard case and / = 1 corresponds to the 
new case. 

As usual we shall make use of a compressibility equation: 

p = p{p) or p = p{p) (16) 

depending on convenience, to describe the physics of the source. 
The boundary conditions at inflnity will be: 

lim S = 1, lim r^E' = (17) 
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where S = A, B or C; and the regularity conditions at the centre of symmetry 
of the configuration will be: 

E', = 0, E'^<oo (18) 

Finally, in those cases where the model has a sharp boundary between an 
interior with p > and vacuum we shall require the continuity of both Xi 
and X', so that the space-time metric will be of global class. 

The system of units we shall use throughout will be such that: 

c=l, 0=1 (19) 

where c is the speed of light in vacuum and G is Newton's constant. A more 
specific system of units of this class will be chosen in Sect. 6. 



2 Explicit field equations 

Taking into account the coordinate condition (|^) and also its derivative: 

B' 6 

C" = 2— - -^{B - C) (20) 

the field equation 6*00 = pA^ can be written as follows: 

BB' B^ 5B'^ B^ . 
^ = -3^-^ + T^ + ^7^-87rpfi-' 21 
rC r^C 4 B r^C^ 

The field equation 5*22 = —ptBCr^ (or 5*33 = —ptBCr^ sin^ 9) can be written 
using the preceding Eq. (O) and the second Eq. (^) as: 



. ,1 , ,A IBA/B \ lAB'^ lA'B' 
A" = ^ -AB'-BA'] + -—- --1] + -—— + 



rC \2 J 2r'^C\C J 8 B^ 2 B 

+47rpB^A + 87TpAB{f/C + B{l-f)) (22) 

The remaining equation we have to take care of is Su = —prB'^, or taking 
into account (|5|): 



B IB'^ B^ A'B' B' BA' 



r'^C r2C2 BA rC rCA 

On the other hand from the conservation equation: 



^ - 2^ = -Snpif + (1 - f)B') (23) 



V^^f = (x^ = r) (24) 
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or: 



s[' + T%{sl - 5°) + vUsl - Si) + vUsl - si) = 

and Eqs. ( [T3| ) we derive the equation: 



(25) 



P' 



p(fB' + l-f)+p +f 



B' 
I 

' B 




P (26) 



Taking into account the regularity conditions (0) at the centre of sym- 
metry of the configuration in the preceding equations it is easy to see that 
they imply: 



Bq = Co, Pto = Pro, Si 
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-BoPto 



(27) 



—PrB^ is satisfied at the ori- 
it follows that it is satisfied 



Therefore the remaining field equation Su = 
gin r = and from from (|21|), ( ^2] ) and ( [2^^ 
everywhere. 

What follows is the summary of this section and the preceding one: The 
models we are considering will be fully described by the field variables A, B, C 
and the source variables p,p; the latter being related by a compressibility 
equation of either type (0). This complete set of variables is constrained 
to satisfy the system of differential equations (|), (0), (P^ ) and (PUI). Ap- 
propriate initial conditions will be Aq, Bq = Cq and po > (or po > 0).The 
initial values of Aq and Bq have to be chosen such that the asymptotic con- 
ditions ( p!7D are satisfied. The boundary of the source will be defined by the 
first zero r = R of the pressure p, and beyond the vacuum field equations 
will be required. The continuity of A, B and C and its first derivatives is 
automatically implemented across the boundary of the source. 



3 Physical and geometrical interpretations 

This will be the more difficult section of this paper, although it does not 
contain any calculation, because it deals about the meaning of words of 
common use. 

When we look at Eqs. (^ as equations to be solved we all refer to the 
r-h-s as the source term. But this is not quite correct because the energy- 
momentum tensor depends on the coefficients of the unknown metric. In 
the case we are considering in this paper the real source variables are the so 
called density p and pressure p related by a compressibihty equation (|TBp. 
The meaning of these three ingredients deserve to be examined with some 
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detail. A density by definition is a mass per unit volume, and a pressure is a 
force per unit surface. Therefore to be clear about them we must tell of what 
mass are we talking about and to what geometry of space are we referring 
when using the words volume and surface. 

As we all know the concept of mass is tricky because it comes in three 
flavours: inertial mass, passive gravitational mass and active gravitational 
mass. Newtonian theory assumes the proportionality of the three masses 
and General relativity assumes the proportionality of inertial and passive 
gravitational masses of test bodies. Beyond that we have a few decisions to 
be taken. 

To decide what geometry of space to use is also a tricky problem in 
relativity theory because we have to decide whether this geometry has to be 
known before we solve the field equations or will be known only after they 
have been solved, in which case the meaning of p, p and the compressibility 
Eq. (|T6|) will also be known only after the problem has been solved. 

Of the three types of mass, only the meaning of active gravitational mass 
was settled very early by Tolman |Q identifying it with the Newtonian mass 
at infinity and proving that it can be calculated as the following integral over 
the source: 

POO 

Ma = An pf"^ df 
Jo 

where we remind that f is the curvature radial coordinate 
formula can be written equivalently as: 

roc 

Ma = ATT pFiry"^ dr (29) 
Jo 

where: 

F(r) = ^y/BC{B'Cr + 25^) 

Bonnor ^ proposed to eliminate the ambiguities 
the passive gravitational mass as: 

/■oo 

M{B)=Ati {p + p)Bf^dr (31) 
Jo 

This means that p + p is interpreted as a density of inertial, or passive, 
gravitational mass and that the metric that gives a meaning to the word 
volume is (0). But this metric is known only once the problem has been 
solved and then it depends on the point of the body which is considered. 
This deprives the meaning of the variables p and p and the compressibility 



(28) 

Using (|^) this 



(30) 

that remain defining 
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equation (|1^) of any a priori significance. Other difficulties that arise from 
this definition are discussed in Bonnor's paper. 

We do not beheve that the consideration of a single body at rest, as we 
have been doing here, can say anything about its passive or inertial gravita- 
tional mass because this would require to know how it reacts to the presence 
of another comparable body. On the other hand we believe that we should 
be able to define its proper mass Mp if we want to know what is the binding 
energy Ei, of any given configuration as defined by: 

E, = Ma- Mp (32) 

More precisely, our point of view, along the lines of a long enduring effort 
to understand the concept of rigidity and establishing a theory of frames 
of reference in special and general relativity, consists in accepting the usual 
generalization of Schwarzschild's "substantial mass" f\ as proper mass: 

Mp = 47i pB^Cr^ dr, Jg = B^Cr^ sin 9 (33) 
Jo ^ 

and defining at the same time a proper mass density pp: 



Pp = pB^C (34) 



such that Mp could be written: 



poo 

p = An / p„r^ dr, = sin 9 (35) 

Jo 

This means then defining density of proper gravitational mass and in- 

terpreting the words volume and surface in the sense of the universal euclidian 
geometry ( |10|) related to the quotient metric (§) by a principal transforma- 
tion (|^). This guarantees the independence of the meanings of Pp, p and the 
compressibility equation (0) independently of the location of the element of 
the fluid in the object and independently of the solution of the field equa- 
tions that one is considering. This guarantees also that Mp can be identified 
with an appropriate number of identical samples of a fluid as weighted with 
a balance at the "shop store" before being assembled into the body. 

The interpretation we have just given of the metric (|10D implies also as 
a corollary that, as above-mentioned, a fluid with isotropic pressure should 
be described by an energy- momentum tensor as written in ([131) . And that 
for the compressibility equation to have a well defined a priori meaning it 
should be given as a relationship between pp and p: 
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An english translation of Schwarzschild's paper has been provided by S. Antoci in 



arXiv:physics/9912033 
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Pp = Ppip) or p = p{pp) 



(36) 



4 Localized energy density 

Following suit to the ideas of the preceding section we exhibit the quantity Eh 
as an integral extended over all space of an energy density function depending 
only on r, B, C and B', C. From (Egj) and (M) it follows that: 



where: 



;>oo 

Eh = 471 pB^Yr'^ dr (37) 
Jo 



while p5 can be obtained from (PT|) as: 



with: 



PB' = - B") (39) 



B{-3rB'-B + ByC) 5 5'^ 



Therefore we have: 



Eh = - XYr'^dr - - / B"Yr^ dr (41) 



2 Jo 2 Jo 

Integrating by parts the second integral we obtain: 



Eh = A'Ki ^(r)r2rfr-47rjim J^fiVM^^ + l- J^l (42) 







1 — >oo 



where: 



with: 



a{r) = ^{XY + Z) (43) 



B' 

Z = ^(~5r'^C^B'^ + 2rBCUVBC + C)B' - 8B^(2VbC-B + C)) 

%rB^^/BC \ J \ II 

(44) 
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depends only on r, B, C and B'. From the asymptotic conditions ( p!7D it 
follows that : 

lim B/C =1 and lim r^B'^ = (45) 

r— >oo r— >oo 

and therefore, the limit in being zero, the final result is: 

Eb = 4:7r / a{r)r^ dr (46) 

5 Linear approximation 

We consider here the linear approximation of the models that we have been 
considering, to take a closer look to two of the de concepts that we have 
implemented in Sect. 3. Namely: the mass defect, or binding energy Eh and 
the proper mass density Pp. 

We assume that A, B and C can be written as: 



A = l + A,, B = l + B,, C = l + Ci (47) 

where Ai, Bi and Ci are small quantities, of order e say. We assume also 
that p is also of order e and that p is of order and can be ignored, as well 
as any other quantity of the same order or smaller, in the field equations. 

The coordinate condition (|^) and the field equations (0), ( [2^ ) become 
then: 



C[ = -{Bi-C,) (48) 
r 

B'l = -3^ + l(5i-Ci)-87rp (49) 

A'l = ^{B[-2A[) + ^^{B,-C^)+Anp (50) 

As our purpose is purely illustrative here we consider below the simplest 
case where the source is a spherical body of finite radius R and constant p. 
The interior solution satisfying the regularity conditions (|18D at the centre 
is: 



2 16 8 

Ai = -npr"^ + ao, Bi = -— vrpr^ + bo Ci = -—npr'^ + bo (51) 
6 15 15 

where and bo are two allowed constants of integration; the exterior solution 
satisfying the asymptotic conditions (|T^ at infinity is: 
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where ai and 63 are two new constants of integration. Demanding the con- 
tinuity of Ai, Bi and Ci and their first derivatives across the border r = R 
fixes qq and all the 6's as follows: 

ao = -27rpR\ 60 = ^vrpi?^ 61 = ^7^pi?^ = ^pR' (53) 

3 3 15 

From the preceding results we can calculate the leading approximation, 

which is of order e^, of the localized energy density (^31). For r < R the result 

is using an arbitrary system of units: 

3 G^MV 

where at this approximation M is either Ma or Mp. For r > R the result is: 

3 G2M2(5r4 + 6i?V-3i?^) 

a = 55 

160 TrrS ^ ' 

The binding energy can be calculated using (|i6|), ( |5^ ) and (|55D, or (p2D, 
(|29|) and (^) at the appropriate approximation. The result is, using arbitrary 
units, the familiar Newtonian amount: 

a result that can be obtained using a variety of other approaches 



6 Massive Neutron Cores 

Any particular model will be characterized by an equation of state and the 
value of its central density, or central pressure, or both in the important 
case in which one assumes that the density is constant. Taking into account 
the regularity conditions the initial conditions of the gravitational potentials 
A, B and C have to be chosen such that: 

Bo = Co, A', = B', = Ci (57) 

And taking into account the asymptotic conditions the values of Aq and Bq 
have to be chosen such that: 

■^See for instance reference H 
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lim A= lim 5 = 1 (58) 

r— »oo r— »oo 

the condition: 

lim C = 1 (59) 

1 — ►oo 

as, well as the remaining asymptotic conditions, being then automatically 
satisfied because the solution behaves as the exterior Schwarzschild one at 
infinity. 

The numerical integration of the system of equations (^, (|2TD, (|2^ ) and 
(|26|) where: 

p = p,{f{B'C)-' + 1 - /) (60) 

and Pp is a known function of p, is a trial and error procedure. Arbitrary 
values of Aq, Bq and po > have to be chosen; the integration has to proceed 
until p = 0', then the equation of state has to be abandoned and p = p = has 
to be required; the integration has then to proceed to sufficiently large values 
of r to check the asymptotic conditions (|58|). If the check is not satisfactory 
the whole process has to be started again with new values of Aq, Bq and 
Po > 0. 

As an important example we consider the equation of state of a degen- 
erate neutron gas as it suits to a model of massive neutron cores. This was 
considered in a famous paper by Oppenheimer and Volkoff from the standard 
point of view which consists in putting / = 0. Our point of view consist in 
using instead the value / = 1. 

The equation of state can be written in parametric form, including both 
points of view, as: 



p = K{smhu-u)if{B^Cy^ + 1- f) (61) 

1 1 
p = -/^(sinhu — 8 sinh -M + 3u) (62) 



where,using arbitrary units: 



4 5 

m c 



m being the mass of a neutron. 

We recall below the values of Ma obtained in []T| for several values of the 
inial values of Uq, and include the values of Mp as calculated from (pSj): 
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Uo 


Ma 


Mp 


Eb 


1 


0.033 


0.033 


-0.000 


2 


0.066 


0.071 


-0.005 


3 


0.077 


0.088 


-0.011 


4 


0.070 


0.085 


-0.015 


5 


0.060 


0.074 


-0.014 


6 


0.049 


0.062 


-0.013 


7 


0.042 


0.054 


-0.012 



The system of units that has been used is that satisfying (JT^) completed 
with the supplementary condition: 

K = l/ATr (64) 

which is the choice made in |]l| to define a unit of mass. The most notori- 
ous result is existence of a maximum mass corresponding approximately to 
uq = 3 whose value is Ma = 0.078 which corresponds to Ma = 0.71 Mq. 
Oppenheimer and Volkoff concluded also from a very crude non relativistic 
argument that above = 3 the equilibrium configurations were not stable. 
For / = 1 the results that we have obtained are: 



Uq 


Ma 


Mp 


Eb 


1 


0.04 


0.04 


-0.00 


2 


0.12 


0.13 


-0.01 


3 


0.24 


0.28 


-0.04 


4 


0.36 


0.47 


-0.11 


5 


0.41 


0.56 


-0.15 


6 


0.36 


0.51 


-0.15 


7 


0.33 


0.48 


-0.15 



Here also we obtain that there is a maximum mass Ma = 0.41 which 
corresponds to Ma = 3.8 Mq but it is substantially larger than the value 
in as well as the mass of some models with anisotropic pressure consid- 
ered by Corchero in . It is even somewhat larger than the limit value, 
Ma = 3.2 Mq obtained by Rhoades and Ruffini in 0] from very general con- 
siderations complying with the conventional point of view. 
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